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Abstract 

We built graphene nanoflakes doped or not with C atoms in the sp 3 hybridization or with Si 
atoms. These nanoflakes are isolated, i.e. are not connected to any object (substrate or junction). 
We used a modified tight binding method to compute the tt and a density of states. The nanoflakes 
are semiconducting (due to the armchair geometry of their boundaries) when their are pure but 
the become conducting when doped because doping removes the degeneracy of the density of states 
levels. Moreover, we showed that the tt Fermi level and the Fermi level of both ir and a electrons 
are not superimposed for small isolated nanoflakes. 

PACS numbers: 73.22.Pr;31.15.bu; 31.15.ae 
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I. INTRODUCTION 



The electronic structure of nanoscale graphene has been the subject of intensive research 
during the last years because of fundamental scientific interest in nanomaterials. With the 
developments in preparation and synthesis techniques carbon-based nanostructures have 
emerged as one of the most promising materials for non silicon electronics. 

Graphene nanoflakes are either semiconductors or conductors aiming to the geometry 
of their boundaries Several studies have been published on graphene nanoflakes doped 
with impurities (2 — 4 1 . We decided here to compute the it and a electrons density of states of 
graphene nanoflakes doped with a small percentage of either Si atoms or with sp 3 hybridized 
C atoms. 

We show here numerical results on the conductivity of graphene nanoflakes and on the 
density of states (mainly the Fermi level) of impurities-doped graphene nanoflakes. 

Our numerical method allows us to take into account the tt as well as the a for the 
calculation of the electronic density of states. And moreover, this method allows one to 
separate the tt and the a density of states. 



II. TIGHT-BINDING MODEL WITH sp" HYBRIDIZATION 

Let us remind that this is a one electron model, each electron moves in a potential V(r) 
which represents both the nuclei attraction and the repulsion of the other electrons, a and 
tt electrons are separately treated. 

The Hamiltonian of the a bonds may be written as: 

H a = E m ^2\iJ >< iJ\ + A a \iJ><iJ'\+Pa \U >< \i'J\ (1) 

where the molecular orbital a is given by: 

\^>=J2 a iJ iJ > ( 2 ) 

i,J 

with \iJ > the hybrid sp u orbital [y = 1, 2, 3) which points from site i along the bond J; 
the energy origin is taken at the vacuum level, i and i' are first neighbours with E m the 
average energy: E m = (E s + vE p )/{y + 1) E s and E p are the atomic energy levels; (3 a is 
the usual hopping or resonance integral in the tight binding theory (interaction between 
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nearest neighbour atoms along the bond); A CT is a promotion integral (transfer between 
hybrid orbitals on the same site): A a = (E s — E p )/{y + 1). 

For an infinite three dimensional crystal (bulk), the gap between valence and conduction 
band (forbidden band) is g — \ — 2/3 a + (v+ 1) A CT | (for VIb elements) and (3 a may be derived 
from the values of g |5j. 

The Hamiltonian for the ir bonds is: 

H w = E p ^2 N >< A +P*^2\i >< *'| (3) 

with \i > the ir orbitals centered on atoms i and (3 n is the hopping integral for 7r levels. 

The (3 n value for C was chosen in order to get the correct positions for C2 energy levels 
in comparison to the results of the Verhaegen's ab initio calculations jfj]. 

We need only three parameters:/?^, (3 n and A a for the homonuclear model which represent 
in fact the average potential V(r) and which take into account the nuclear attraction and 
the dielectronic attraction [5[ . The values of these three parameters are given in table I. 

III. RESULTS 

We computed the density of states of pure and doped graphene nanoflakes. The nanoflakes 
have boundaries represented in figure 1. We studied nanoflakes containing 480, 720 and 960 
atoms: we obtained their geometry by reproducing figure 1 in the y direction. 

In figure 2, one may see the density of states of nanoflakes of graphene containing 480 
atoms. The nanoflakes are isolated, i.e. they are not connected to any other conducting 
material or connected to any substrate. The 3 graphs in figure 2 representing the DOS 
(density of states) correspond to: (a) pure graphene nanoflake (b) 5% of Si atoms randomly 
introduced in the graphene nanoflake and (c) 5% of sp 3 hybridized carbon atoms randomly 
distributed within the nanoflake. , 10%, 15% and 20%. The boundaries of the nanoflakes 
correspond to the armchair geometry. 

Figure 3 corresponds to the same nanoflakes (i.e. pure, with 5% of Si atoms and with 
5% of sp 3 hybridized carbon atoms) but in this case for 720 atoms. 

Figure 4 corresponds also to 3 DOS with the same nanoflakes but this time containing 
960 atoms. 
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It is possible to introduce impurities like Si atoms. Carbon and Silicon belong to the 
same column IVb of the periodic table. However, their structures in the solid states are dif- 
ferent. Carbon indeed can be found in several allotropic varieties such as diamond, graphite, 
graphene whereas silicon is always found in three dimensional systems. 

In mixed silicon-carbon Si n C p nanoclusters, it has been previously pointed out that 
experimental results show a progressive transition from a carbon type behaviour of Si to 
a silicon type behaviour as the silicon richness is growing in the clusters jg|. This may 
be transferred in nanoflakes to the percentage of Si: here we deal with small percentages 
of Si surrounded by 95% of C at least. From references Q] we may say that due to the 
richness in carbon atoms in the studied nanoflakes of graphene, the Si atoms adopt the 
carbon behaviour i.e. are in the sp 2 hybridization. 

IV. DISCUSSION 

It is worth noticing that for small and finite size nanoflakes, the DOS consists of many 
discrete levels and instead of Fermi level, one should speak of highest occupied molecular 
orbital (HOMO) and lowest unoccupied molecular orbital (LUMO). 

However, we retain these terms in a loose sense because the main result when comparing 
the three figures of DOS for the nanoflakes of graphene of different sizes, we observe (figures 
2,3 and 4) that the Fermi (or HOMO) level of the it electrons DOS correspond to the DOS 
in the k = direction in the Brillouin zone (see for that reference 9j. The fact is that 
for 480 atoms and 720 atoms of pure graphene (figures 2(a) and 3(a)) the DOS shows that 
there is a gap between the ir Fermi level and higher ir energy levels: this indicates that the 
nanoflake of pure C atoms is semiconducting. This is in good agreement with the literature 
showing that with armchair boundaries (which is the case in our nanoflakes) semiconducting 
is generally obtained [lj. 

By doping the nanoflake with either Si atoms or by sp 3 carbon atoms, one destroys the 
semiconducting character. This is due to the fact that impurities remove the degeneracy of 
the energy levels. 

What is also interesting is that the 7r Fermi level and the whole Fermi level are not 
superimposed for nanoflakes containing a small number of atoms (480 and 720 atoms: fig. 
2 and 3). The superimposition only occurs for a larger number of atoms (here 960 atoms: 
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fig.4). Generally in the literature, one studies the ir DOS as the major numerical methods 
used by authors is the tight binding one taking only account of n electrons. But here, we 
deal with a and n atoms, and when calculating the Fermi the n and the whole Fermi level, 
these two do not superimpose, at least for a small number of atoms. But a electrons do not 
enter in the conducting or semiconducting character of the sample. 



V. CONCLUSION 



Let us summarize the results: doping graphene nanoflakes with armchair boundaries 
removes the degeneracy of the energy levels in the DOS and therefore removes the semi- 
conducting character of graphene nanoflakes with this kind of boundaries; the n Fermi level 
and joint n and a Fermi level do not superimpose for small graphene nanoflakes: but as the 
a electrons do not enter conduction this has no consequence on the electronic features of 
small graphene nanoflakes. 
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atom /parameter 


E a 


E n 


Ax 


Ar 


C 


19.45 


10.74 


7.03 


3.07 


Si 


14.96 


7.75 


4.17 


0.8 



TABLE I: Parameters for the tight binding calculations in eV 



FIG. 1: Example of graphene nanoflake with 240 C atoms. Our graphene nanoflakes are obtained 
by adding the same structure in the y direction 
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FIG. 2: Electronic density of states of a graphene nanoflake containing 480 atoms, (a) pure c; (b) 
doped with 5% of Si atoms; (c) doped with %5 of C atoms in the sp 3 hybridization 
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FIG. 3: Electronic density of states of a graphene nanoflake containing 720 atoms, (a) pure c; (b) 
doped with 5% of Si atoms; (c) doped with %5 of C atoms in the sp 3 hybridization 
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FIG. 4: Electronic density of states of a graphene nanoflake containing 960 atoms, (a) pure c; (b) 
doped with 5% of Si atoms; (c) doped with %5 of C atoms in the sp 3 hybridization 
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